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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. We present the results of a detailed spectroscopic analysis of 20 years of observations of 
AG Carinae using the radiative transfer code CMFGEN. Among the conclusions of this work, we highlight the 
importance of including time-dependent effects in the analysis of the full S Dor cycle. We obtained that the 
mass-loss rate is approximately constant during the cool-phases, implying that the S Dor-type eruptions begin 
well earlier than the maximum seen in the visual lightcurve. We also determined that the S Dor cycles are 
ultimately a consequence of an increase/decrease of the hydrostatic radius in combination with the formation 
of a pseudo-photosphere. 

ABSTRACT 

We present the results of a detailed spectroscopic analysis of 20 years of observations of AG Garinae using 
the radiative transfer code CMFGEN. Among the conclusions of this work, we highlight the importance of 
including time-dependent effects in the analysis of the full S Dor cycle. We obtained that the mass-loss rate 
is approximately constant during the cool-phases, implying that the S Dor- type eruptions begin well earlier 
than the maximum seen in the visual lightcurve. We also determined that the S Dor cycles are ultimately a 
consequence of an increase/decrease of the hydrostatic radius in combination with the formation of a pseudo- 
photosphere. 
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1. AG CAR IN CONTEXT 

The massive star AG Carinae has been the most 
variable object among the Galactic luminous blue 
variables during the last 20 years. Dramatic pho- 
tometric changes of 3 magnitudes in the V band 
have been recorded, while the stellar spectra have 
indicated effective temperatures ranging from about 
9kK to 25 kK. AG Gar is a rare member of the LBV 
class that shows strong S Dor-type variations and is 
surrounded by a bipolar nebula. 

This work is a result of a monitoring campaign 
spanning 20 years of observations of AG Gar to ob- 
tain its stellar and wind parameters for each epoch. 
The first results were discussed in Groh et al. (2006), 
highlighting the first observational evidence of fast 
rotation in LBVs. The evolution of both stellar (ef- 
fective temperature, radius, luminosity) and wind 
parameters (mass-loss rate and wind terminal ve- 
locity) along the S Dor cycle, and the correlations 
among them, provides significant insights about the 
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physical processes going on during this short and un- 
stable evolutionary stage. 



2. SPECTROSCOPIC MODELING 

The physical parameters of AG Car were ob- 
tained through spectroscopic modeling using the ra- 
diative transfer code CMFGEN (Hillier & Miller 
1998), which allows a simultaneous solution of the 
radiative transfer and rate equations. Spherical ge- 
ometry and steady state are assumed to solve those 
equations, while clumping is included using a depth- 
variable filling factor. Fully-blanketed models were 
calculated for AG Car, taking into account the pres- 
ence of over than 50,000 transitions of H, He, C, N, 
O, Si, P, Na, Mg, Fe, Mn, Co and Ni. 

For the first time in the analysis of LBV winds, 
time-dependent effects were included consistently in 
the velocity law and density structure of the wind 
assuming that the mass-loss rate and wind terminal 
velocity change linearly as a function of time (J. H. 
Groh et al. 2007, in preparation). Significant differ- 
ences can be seen in the fits obtained with a time- 
dependent model in comparison with a steady-state 
model, especially in the ultraviolet region (Fig.[T]). 



1 



2 



GROH, DAMINELI, & HILLIER 



TABLE 1 

STELLAR PARAMETERS DURING MINIMUM 

Epoch Toff M Woo 
{Rq) (K) (Mpyr-i) (kms-^) 

1985-1990 62 23,800 1.5 x 10"^ 300 
Apr 2001 69 20,500 2.5 x 10"^ 105 



phases, significant differences are derived when time- 
dependent effects are taken into account. In spe- 
cial, M is actually approximately constant during 
the cool phase (1991-1999), in comparison with a 
strong increase in M around maximum when steady- 
state models are used. Thus, the results obtained 
from time-dependent modeling imply that the erup- 
tions during the S Dor cycles begin well earlier than 
the maximum seen in the visual and near-infrared 
lightcurve. In addition, it is misleading the idea of 
an increase in M as a consequence of the maximum 
- actually, the maximum in the lightcurve is a con- 
sequence of the increase in M which occurred 2 or 3 
years ago, in the beginning of the cool phase. 

3.3. Changes in the photospheric and hydrostatic 
radius 

The evolution of the photospheric and hydro- 
static radius of AG Car during the S Dor cycle is 
presented in the panel b of Figure [H It is apparent 
to conclude that the S Dor-type variability is driven 
mainly by changes in the hydrostatic radius of the 
star. Changes in the photospheric radius do occur 
as a result of the dense wind and variable mass-loss 
rate, but this change solely is not sufficient to explain 
the strong temperature changes seen from minimum 
to maximum. 
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Fig. 2. Evolution of the mass-loss rate (panel a) and 
of the hydrostatic and photospheric radius (panel b) of 
AG Car during the S Dor cycle. 
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3. RESULTS AND CONCLUSIONS 

3.1. Comparison between the minimum epochs 

The stellar and wind parameters are different 
during consecutive minima as derived from the anal- 
ysis of spectra taken in 1985-1990 and 2001 (Ta- 
ble [1]). During those epochs, a significant difference 
of ~ 3000 K in Toff is noticed. In addition, a remark- 
able difference is seen in M (which is 2 times higher 
in 2001 than in 1985-1990) and Voo (3 times lower 
in 2001). We interpret the changes in M and Vao as 
a consequence of the bistability of line-driven winds 
(Lamcrs et al. 1995). Indeed, there is a change in the 
iron ionization structure around the sonic point, as 
predicted by Vink & de Koter (2002). The dominant 
iron ionization stage below the sonic point (which 
is around 18 kms~^) changes from Fe^+ in 1990 to 
Fe2+ in 2001. 
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Fig. 1. Comparison between the ultraviolet spectrum 
of AG Car taken in 1991 October (grey) and a steady- 
state model (panel a, dashed) and a time-dependent 
model (panel 6, dashed) . The better quality of the time- 
dependent fits is remarkable. 



3.2. The mass-loss rate evolution 

Figure [5] (panel a) shows the mass-loss rate evo- 
lution along the S Dor cycle obtained from steady- 
state and time-dependent models. As a result of 
the large flow timescale 5 years) during the cool 



